Universal Markovian reduction of Brownian particle dynamics 
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Non-Markovian processes can often be turned Markovian by enlarging the set of variables. Here 
we show, by an explicit construction, how this can be done for the dynamics of a Brownian particle 
obeying the generalized Langevin equation. Given an arbitrary bath spectral density Jo, we intro- 
duce an orthogonal transformation of the bath variables into effective modes, leading stepwise to a 
semi-infinite chain with nearest-neighbor interactions. The transformation is uniquely determined 
by Jo and defines a sequence { J n }neN of residual spectral densities describing the interaction of the 
terminal chain mode, at each step, with the remaining bath. We derive a simple, one-term recur- 
rence relation for this sequence, and show that its limit is the quasi-Ohmic expression provided by 
the Rubin model of dissipation. Numerical calculations show that, irrespective of the details of Jo, 
convergence is fast enough to be useful in practice for an effective Markovian reduction of quantum 
dissipative dynamics. 



Introduction. As is well known, the study of open sys- 
tems, in both the classical and quantum case, is a sub- 
ject of major interest in physics, chemistry, and various 
other disciplines. In many applications and fundamental 
experiments, one is faced with the reduced dynamics of 
a relatively simple subsystem which can be manipulated 
and measured, while the environment is only partially 
under control. A thorough understanding of the ensuing 
dynamics has been obtained for the Markovian case, in 
which feedback from the environment to the system can 
be neglected, and general analytical results are available 
together with efficient numerical algorithms [H-Q- The 
situation is much more involved in the non-Markovian 
regime, which typically arises due to strong coupling 
and similar time scales of system and bath evolution. 
In this case general strategies are still available 0, 01 , 
but differently from the Markovian case they typically 
lack simple results of general validity, to be expressed 
in terms of the phenomenologically relevant quantities 
and leading to manageable numerical tasks. A bridge 
between the two situations can be built relying on a 
suitable embedding of a non-Markovian dynamics in a 
Markovian one, as recently addressed in [4j-[6|| . Indeed, 
while it is common wisdom that a non-Markovian pro- 
cess can be embedded in a Markovian one by a suitable 
enlargement of the number of relevant variables already 
at classical level Q , there is no universal recipe for how 
this can be done and which class of non-Markovian pro- 
cesses can be reached. 

In the present Letter we demonstrate how such a 
Markovian reduction can be achieved for the ubiqui- 
tous model of quantum dissipation provided by a Brow- 



nian particle, or a two-level system, linearly coupled to 
a bath of harmonic oscillators characterized by an arbi- 
trary spectral density (SD) Q. The procedure is phys- 
ically transparent, in that it focuses exclusively on the 
SD, and all relevant quantities can be constructed in 
terms of the SD. As will be shown below, the system 
dynamics can equivalently be described including, be- 
sides the Brownian particle degree of freedom, a set of 
effective environmental modes coupled in a linear-chain 
fashion. The terminal mode of the chain couples to a 
residual bath and undergoes a Brownian-like dynamics 
which rapidly approaches a Markovian behavior over 
the whole interval of relevant frequencies as the length 
of the chain increases. The model as such is cl osely re- 
lated to Mori's theory Q and its generalizations (ToT 11 1 . 



While previous work by two of us 0, @| has focused on 
the implications of a Markovian truncation of such ef- 
fective mode chains, the present analysis proves the con- 
vergence towards Ohmic behavior, and thus the general 
validity of the procedure. The question of how to cor- 
rectly set the initial state of the chain will be detailed 
in a forthcoming paper. 

Effective-mode transformation. We start by consider- 
ing the Caldeira-Leggett Hamiltonian, here written in 
mass- weighted bath coordinates x^, 
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which is known to lead, in the continuum limit, to a gen- 
eralized Langevin dynamics for the system described by 
the s degree of freedom. The reduced system dynamics 
is entirely determined by the SD of the environmental 
coupling Jo(w) which, for the microscopic model above, 
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reads as Q 
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In general, Jq(uj) is a real, odd parity function defined 
by the real part of the frequency-dependent memory 
kernel [ItJ 7(w) entering the generalized Langevin equa- 
tion (GLE), namely Jo(u>) = mu>Bjej(u;) and Jo(u>) > 
for u) > 0. It fully determines 7(0;) by virtue of the 
Kramers-Kronig relations, as well as the correlation 
function of the GLE random force by virtue of the 
fluctuation-dissipation theorem. In the following we as- 
sume, as a typical situation, that Jq(lo) is strictly pos- 
itive and continuous in an interval (0,Wij) - where lor 
is a high-frequency cutoff - and zero otherwise on the 
positive real axis; other interesting cases will be briefly 
considered below. 

Given a GLE and its relevant SD Jo, Eq.((2]) allows to 
define a microscopic model for the dissipative dynamics 
of the s degree of freedom, e.g. by introducing a bath 
of harmonic oscillators with evenly spaced frequencies 
Uk = kAuj (k = 1, . . . , N) and setting the coupling co- 
efficients of Eq.([T]) as 
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The system-bath interaction term in Eq.([T|), H mt = 
— ^2^ =1 CkXkS = —DqXis, naturally introduces an ef- 
fective mode X\ = ~^2k=i c k x k/Do where Do is a nor- 
malization constant which in the continuum limit reads 
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if Jo°° dcoJo(io)oj. This defines the first 
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column of an otherwise arbitrary, orthogonal matrix T 
transforming the bath coordinates x l = (x\, . . . , xjy) 
into X* = {Xi,X' 2 ,...,X' N ), X = T*x. The trans- 
formation can be fixed by requiring that the "residual" 
bath of coordinates X' 2 , X' 3 , . . . , X' N is in normal form, 
i.e. that (T t u) 2 T) lj = SijQ 2 holds for all i,j > 1 with 



= 1,...,N) . Note that the fre- 



quency fij of the effective mode introduced in this way is 
solely determined by the SD J (w), O. 2 = (T*u; 2 T)ii = 
T,k=i^k c l/ D « ^*fi° du Mu)u a - The couplings 
C k = -(T*w 2 T)i ifc (k = 2,...,N) between the nor- 
mal modes of the residual bath and X\ allow one to 
re- write Eq.JTJ) as a Caldeira-Leggett-like Hamiltonian 
for the Xx degree of freedom, thereby defining the SD 

Ji(u) = § E* =2 j£<*(w ~ ^ fe ) " felt " by the effective 
mode Xi, which is the only bath mode directly cou- 
pled to the s degree of freedom. Clearly, in the con- 
tinuum limit, the procedure can be indefinitely iter- 
ated and used to define a sequence of effective modes 
Xi, X2, ■ ■ ■ , Xm , ■ ■ ■ coupled in a linear-chain fashion 
and a corresponding sequence of SD Ji, J 2, ■ ■ ■ , Jm, ■ ■ ■ 
characterizing the residual bath "felt" by each mode, 
see Fig|TJ In other words, there exists an orthogonal co- 
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Figure 1: (Color online) Schematic for the linear-chain trans- 
formation described in the text. X„ are the effective modes, 
On their frequencies and D n the couplings between adja- 
cent modes X n -i and X n (X = s). The SD Jm (M > 0) 
describe the interaction with the residual bath. 



ordinate transformation which converts the continuum 
version of Eq. ([I]) into the form 
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where, for n 0, 
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(and similarly for P n +\ in terms of p(uj)); AV(s) = 
5Q 2 s 2 /2 is a counter term involving the renormaliza- 
tion frequency <50q = (2/w) J Q dui Jq(u))/oj. As we show 
below, though different canonical transformations to a 
linear chain be devised (see e.g. Ref. [T2I ) . the one sug- 
gested above, with the coupling in pure coordinate form, 
allows one to write an explicit expression for {J n }neN 
without knowing the eigenfrequencies of the residual 
bath at each step. 

Sequence of SD. As observed by Leggett 
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the SD acting on the system degree of freedom can 
be obtained from the analytically continued, Fourier- 
transformed classical (or Heisenberg) equations of mo- 
tion as the imaginary part of a propagator, 
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This procedure has recently been used by two of us to 
obtain a continued-fraction expression for the SD gen- 
erated by a linear chain with Markovian closure, which 
in turn formed the basis for approximating a given SD 
Employing a similar strategy, we now focus on the 
properties of the residual spectral densities Jm closing 
the chain after M effective modes have been extracted 
as outlined above. For the Hamiltonian of Eq. ([1]) , after 
introducing the first effective mode Xi, we obtain 
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where tik and Ck have been introduced above. In the 
continuum limit, with the help of Eq.Q, the sum in the 
denominator can be replaced by the function [18 1 
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or, equivalently, 
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In this form W\ is given as an integral of its limiting 
imaginary part, Ji(u>) = ImW^ (u>) . In the following, a 
function W\ defined by E q.(pt will be referred to as the 
Cauchy transform of J\ [la ]: it is an analytic function 
in the whole complex plane except for the support of J\ 
on the real axis [19(, which vanishes as z~ 
We define 
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which, analogously to L = — z 2 — W , gives Jq{lo) = 
IhiWq (uj). It follows that J\ can be written in terms 

of Jo as Ji(cj) = ^rr\k where, as we now show, Wq 
is the Cauchy transform of Jo. In order to prove this, 
we notice that according to its definition, EqJBJ), Wo 
is analytic in the upper and lower half planes [20(, and 
vanishes as z~ 2 for \z\ — > 00. Writing Wo(z) as a Cauchy 
integral on a large semicircle in the upper half plane, 
we can add a term ±(lj — z*)" 1 to the integrand, and 
get, from the real and imaginary parts of the resulting 
expression, the desired result. In general, then 
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defines a recurrence relation for the SD J n+ i(u>) felt by 
the n + 1-th effective mode, given the SD J n (oj) of the 
n-th mode. Equivalently[2l|, 
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is a recurrence relation for the Cauchy transforms 
which only requires the first Cauchy transform (Wo) 
as an input and easily provides the sequence J n (ui) = 
TmW+(u). 

Eq.© represents the main result of this Letter. It is 
a simple recurrence relation between the Cauchy trans- 
forms of the SDs which allows us to write the limiting 
condition as 



W(z) = Cl 
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provided Q = lim n Q n and D = lim„ D n exist. The 
physical solution (ImW + ^ 0) provides the SD "closing" 
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Figure 2: (Color online) Left: the SD J defined in Eq.{TT]) 
for cj = 0.04 a.u., do = 0.01 a.u. and 7 = 0.01 a.u.. Middle: 
results of the deconvolution of Jo for the first five modes, 
obtained when setting the high-frequency cutoff ujr to 0.1 
a.u.. The Rubin SD of Eq. (|10p with the same ujr is shown 
as dashed line. Right: Effective modes parameters (O^ and 
D n ) up to n — 10. 



the chain, which has a non-vanishing value for u> 2 G 
[fl 2 — 2D, fl 2 + 2D] only. In other words, the limiting 
SD reads as 

J(ui) = i^/ (uj 2 - uj 2 l )(uj 2 r - uj 2 ) uj l < u < uir (9) 

where u? R = fl 2 + 2D and lo\ = Max{Cl 2 - 2D, 0}. The 
requirement Jq{oj) > for U) > fixes fl 2 = 2D, since 
the condition Cl 2 > 2D (though physically admissible) 
would give rise (using the above recursion procedure 
backwards) to a SD Jo with a low- frequency cutoff ujl = 
Vfl 2 - 2D. Therefore Q, O 2 = 2D, uj 2 r = AD = 20 2 
and Eq.® reduces to the quasi-Ohmic SD provided by 
the Rubin model of dissipation 0, [l6j , 
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This means that provided a sufficient number of effec- 
tive modes is included in the definition of the system, 
the resulting dynamics is Markovian. In practice, as 
we show numerically below this number is rather small, 
since convergence is quite fast even for structured spec- 
tral densities. Notice though that when Jo(uj) has a 
low frequency cutoff lol but is otherwise positive on the 
interval (ujl,^r), Eq.© shows that no Markovian re- 
duction is possible, no matter how many effective modes 
are included in the system. 

In general, a high-frequency cutoff lur can be natu- 
rally associated to the SD Jo(^), determining the spec- 
trum of environmental frequencies relevant for the re- 
duced system's dynamics. This suffices to show that 
ReW + (w) diverges logarithmically for ui — > ojr, unless 
^0 (w.r) is equal to zero. In view of Eq. ([7|) this in turn im- 
plies that Ji(w) -> as w -> lur, and this cutoff is later 
on automatically preserved. Furthermore, starting from 
Eq. ([8]) one can also immediately obtain some interesting 
bounds on the value ReW^"(oj) can take at the extreme 
points of the relevant frequency interval (0,w^), which 
help determine the behavior of the recurrence relation 
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Eq.©. Among these we only note here that ReVF+(w) 
is positive for ui = and negative for to = lur, and then 
by continuity the function will go through zero at some 
intermediate point u). This explains why the procedure 
generally fails to converge for SD with gaps, since if 
H) falls in the gap one has W^(ui) = and, by virtue 
of Eq.©, this introduces an isolated pole in W n +i{z) 
which invalidates the use of Eq.©. Notice, however, 
that even in this case an orthogonal transformation of 
bath variables into linear-chain modes can still be intro- 
duced to define a number of chains of effective modes, 
one for each interval where Jo(uj) > 0. 

Numerical results. The procedure described above 
for the determination of the sequence { J n }neN of effec- 
tive SD can be easily implemented numerically, relying 
on the recurrence relation Eq.(j8]) for the SD Cauchy 
transforms [23], the only necessary input being the ini- 
tial SD Jo and the cutoff ujr. To show the effectiveness 
of the method and rapidity of convergence we consider 
the numerical results for a couple of representative SDs, 
with a frequency cutoff ojr fixed in such a way that 
Jo(w/j) « 0. As a first example we consider 

Jo(u) = §TT~ i — T~ 2' (11) 

(or — loqY + j z uj z 

which is the effective SD felt by a Brownian particle 
coupled to a harmonic oscillator of frequency ujq which 
in turn interacts with an Ohmic bath[li]. This cou- 
pling scheme is evident from the results of Fig |2] where 
Ji, plotted in the middle panel, appears to be Ohmic, 
as can also be checked analytically. The sequence then 
very rapidly becomes indistinguishable from the Rubin 
SD given by Eq.fllUp for the chosen cutoff u)r, as can also 
be seen from the ratio fi^/£) n — > 2. As a second exam- 
ple we consider a highly structured, multipeaked SD as 
plotted in Fig(3l It is clear from the figure that also 
in this case convergence is quite fast, and the limiting 
Rubin SD is obtained after few (say 10-15) iterations. 

Conclusions. We have presented a recursive proce- 
dure to recast the non-Markovian dynamics of a Brow- 
nian particle, interacting with a bath characterized by 
an arbitrary SD, into the Markovian dynamics of an en- 
larged set of variables including effective modes of the 
reservoir coupled to a quasi-Ohmic residual SD. The ap- 
proach provides an explicit analytic relationship among 
successive residual SD, which can be easily evaluated 
numerically starting from an arbitrary (gapless) initial 
SD. These results pave the way for an efficient general 
treatment of quantum dissipation in the presence of ar- 
bitrarily complex environments. 
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